THREE-FREQUENCY RESONANCES 
IN DYNAMICAL SYSTEMS 
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We investigate numerically and experimentally dynamical systems having three interacting frequencies: 
a discrete mapping (a circle map), an exactly solvable model (a system of coupled ordinary differential equa- 
tions), and an experimental device (an electronic oscillator). We compare the hierarchies of three-frequency 
resonances we find in each of these systems. All three show similar qualitative behaviour, suggesting the 
existence of generic features in the parameter-space organization of three-frequency resonances. 



1. Introduction 

A fundamental goal in the modelling of complex be- 
haviour by means of low-dimensional nonlinear de- 
terministic systems is the understanding of the tran- 
sition to chaos from quasiperiodic motion on a T-* 
torus for small j. This transition is well understood 
for the case j = 2, the paradigmatic model being the 
periodically forced nonlinear oscillator and its dis- 
crete version: the circle map. 

In T^, orbits can be hierarchically arranged by 
means of rotation numbers [Arrowsmith & Place, 
1990] 
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periodic orbits — lockings or resonances — being 
then characterized by rational rotation numbers. If 
two rational rotation numbers (p/g, r/s) satisfy \ps - 
rq\ — 1 they are said to be unimodular or Farey adja- 
cents. Between two periodic solutions characterized 
by unimodular rotation numbers there exists a pe- 
riodic solution with a minimal period. This solution 
is given by a rotation number termed the Farey me- 
diant of the two previous ones: m/n = p/q ® r/s = 
{p + r)/{q + s). Moreover this periodic solution is the 
most prominent — the largest — in the open interval 
{p/q, r/s) between the two parents. Continued ap- 
plication of the mediant gives the Farey tree which 
underlies the organization of the locking intervals of 
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Figure 1: (a) Three-frequency devil's staircase in the quasiperiodically forced circle map. The external 
frequencies are A = 1 and A = 12/7. The steps at left and right correspond to the plateaux (12/7)/7 = 
0.2448979 ... and 1/4 = 0.25; the largest plateau between the two is (1 + 12/7)/ll = (19/7)/ll = 0.2467532 . . ., 
in accordance with the prediction of the generalized Farey mediant, (b) Devil's ramp system: the rotation 
number as a function of the external frequency ratio and the intrinsic frequency for a quasiperiodically forced 
circle map representing a doubly kicked nonlinear oscillator. The devil's ramps illustrate the overall ordering 
of three-frequency resonances. 



rational rotation numbers in the space of parameters 
of tw o-frequency systems; the so-called devil's stair- 
case []Gonzalez & Piro, 1983t lAronson et al., 1983; 



1989; Arrowsmith et al., 1993]. These properties of 



Gonzalez & Piro, 1985; [Cvitanovic et al., 1985 ; Hao, 



rotation numbers are closely related to the contin- 
ued fraction expansion which gives the best rational 
approximants to an arbitrary real number. 

In order to study torus breakdown we have 
generalized this approach to the case of the simul- 
taneous approximation o f a pair of real numbers 
[ Cartwright al., 1999b ]. Three-frequency systems 
also possess a structure of resonances, but this time 
more complex, for in addition to the rational rela- 
tions found in two-frequency systems there is now 
a new type of locking: a three-frequency resonance, 
given by the nontrivial solutions of the equation 
nfi + m/2 + qfs = with n, m, and q nonzero inte- 
gers. Three-frequency resonances form a web in the 
parameter space of the frequencies [Baesens et al. , 
1991]. 



We have found that [Cartwright et al., 1999b]: 



• We can define a subharmonic real interval 
(/i/p, /2/g) inside which the hierarchical orga- 
nization of three-frequency resonances is well 
described by a generalized Farey sum. 

• In order to implement the generalized Farey 
sum we must modify the adjacency condition: 
if p/g is a convergent of the forcing frequency 



ratio we define two fractions of real numbers 

fx/m, fy/n as adjacents if lUn - fym\ = \ fiq - 
f2p\ where /i and /2 are the external frequen- 
cies. 

• We can now define the generalized Farey sum 
between fractions of real numbers which satisfy 
this adjacency condition: fr/k = fx/m (S fy/n = 
{f. + fy)/{m + n). 

Here we present numerical and experimental 
studies of three dynamical systems having three in- 
teracting frequencies, in which we can observe three- 
frequency resonances, and compare our observations 
with the above predictions. 



2. Numerical Simulations 

We have investigated a discrete mapping: a 
quasiperiodically forced circle map given by the 
equation 



k 

+ f2„ sin 27: (j) mod 1. 
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(2) 



The quasiperiodic sequence il„ is the time interval 
between successive pulses of a sequence composed of 
the superposition of two periodic subsequences, one 
of period Ti = 1 and the other of period T2 (T2 < 1 
with no loss of generality), multiplied by the value of 
the intrinsic frequency of the oscillator. 



(a) 



(b) 




Figure 2: (a) The frequency response of both oscillators of Eq.(^) as a function of the intrinsic frequency of the 
first oscillator. Zone I is chaotic, with intermittent synchronization, in zone II there is 1/1 synchronization, 
and in zone III, synchronization other than 1/1. The input frequencies have been set to be equivalent to those 
of Fig. |l[ (b) A magnification of the central zone of (a), where the oscillators are mutually synchronized at 1/1. 



In Fig. I^a we can observe that there are intervals 
of constant frequency response, namely the horizon- 
tal segments of the figure. These intervals represent 
stable three-frequency resonances, that is, responses 
that satisfy nfi + m/2 + qfs = with n, m, and q 
nonzero integers. We can also see in Fig. |^a that the 
stability widths of the resonances form a hierarchi- 
cal structure very similar to that of the well known 
devil's staircase in periodically forced nonlinear os- 
cillators. This is the generalized devil's staircase for 
three-frequency systems. 

The other image. Fig. j^, shows the devil's 
ramps. These represent the global hierarchy of 
three-frequency resonances shown as a function of 
the external frequency ratio and the intrinsic fre- 
quency for a critical {k = 1, the map instantaneously 
has critical behaviour) quasiperiodically forced circle 
map. 

We have also performed numerical simulations 
with two parametrically coupled nonlinear oscilla- 
tors, each with an exact analytical solution, forced 
by means of two impulsive periodic forces. The differ- 
ential equation for each osci llator [Gonzalez & Piro, 
1983; [Gonzalez & Piro, 1985P is: 



Ui + {Abui-2a)u, + b^ul-2abui + {ijJar{tY + a^)ui = 

(3) 

for i = 1,2. The coupling and forcing terms woi(i) = 
(Doi + sgnUi{t)sgnUj{t)Ai and fi{t) = VsiJ^n^i^ ^ 
nTEi) preserve the piecewise integrability of the sys- 
tem. 

We have made a power spectrum analysis of the 
output of both oscillators. In Fig. ^a we display the 
most prominent peak in each spectrum versus the 
intrinsic frequency of oscillator 1 for a parameter 
region equivalent to that of Fig. |l|a. In Fig. ^ we 



show a magnification of the zone of Fig. ^ in which 
the two oscillators are synchronized at 1/1. The 
three principal peaks in the Fourier spectrum satisfy 
nfi + m/2 + qfi = with n,m, and q nonzero integers. 
All the other peaks in the spectrum can be expressed 
as linear combinations of this fundamental set. 



3. Experimental Results 

We have constructed an electronic oscillator — Fig. ^ 
— a higher-dimensional version of a phase-locked 
loop (P LL), forced with t wo independent periodic 
forces [ palvo et al., 1999 ). Our circuit consists of 
two coupled voltage-controlled oscillators forced with 
two external forces of frequencies /i and /2. As 
a basic circuit for both oscillators we use a digital 
phase-locked loop integrated circuit, the CD 4046A. 
The outputs of the two phase-locked loops are sent 
to a type 1 phase comparator. The error signal is 
fed back to both voltage-controlled oscillators, and 
passes through an overall adjustable amplifier to 
provide control over the coupling strength; we are 
interested in the weak coupling regime in this work. 
Inverted and direct versions of the error signal are 
sent to oscillators one and two respectively; this in- 
version of the error signal in one of the paths is 
necessary for the stability of the circuit. Feedback 
signals enter the voltage-controlled-oscillator control 
pins through appropriate adder circuits. The adders 
also allow independent coupling with the external 
forces and tuning of the internal frequencies through 
application of adjustable DC levels. 

In Fig. ^ we plot the fundamental frequency /a 
versus the DC offset of one of the forcing signals. As 
in the previous numerical simulations, the general- 




Figure 3: The circuit. The output of the voltage controlled oscillator (VCO) sections of both digital phase-locked 
loop devices are sent to an exclusive OR port for phase comparison. The phase comparator output is returned, 
after low-pass filtering, to both VCO inputs. An amplifier in the feedback path controls the coupling strength 
between the oscillators. Appropriate operational adders on both VCO inputs allow the external forcing of the 
oscillators. Stable phase-locked responses between oscillators also require an additional unit gain inverter 
prior to one of the device inputs. 



ized mediant is the largest plateau between the two 
adjacents given by the subharmonics of the forcing 
frequencies. The circuit of Fig. ^ represents qualita- 
tively a large class of real dynamical systems charac- 
terized by some kind of feedback loop with a trans- 
ference function of the low-pass type, and can be 
useful in clarifying the physical mechanisms that 
underlie three-frequency resonances. In Fig. ^ we 
show the Fourier spectrum of the oscillator response 
at a point inside the stability region for the three- 
frequency resonance. In Fig. ^b, a magnification 
around this frequency value, we can see that har- 
monics of the main resonance at fs = 518.2 Hz lie 
at equal distances S = 27.3 Hz from both external 
frequencies. The peak at fs dominates the spectrum, 
being greater than those corresponding to the forcing 
frequencies. The distance S = sfs appears as a mod- 
ulation of the spectrum, which shows, on this fine 
scale, a series of peaks separated precisely by this 
distance. 

The origin of these frequencies can be interpreted 
in terms of the general structure of the system. Pas- 
sive nonlinearities are able to generate appropriate 
frequency harmonics. These, in turn, through non- 
linear mixing with the driving frequencies, can gen- 
erate the frequency S. Only for the case of the main 
three-frequency resonance are the two Ss of the same 
value, i.e., 6 = Si =62- Thus, further nonlinear mix- 
ing in this case gives only terms of zero frequency 
and harmonics of S, that is, no other frequencies are 
added to the system. Otherwise, for an arbitrary 
frequency, the two Ss are in general different and, 



if irrationally related, generate an infinity of low- 
frequency components that pass through the low- 
pass filter and are fed back to the system, destabiliz- 
ing the response for this particular frequency value. 
This physical mechanism can also explain the order- 
ing of the three-frequency resonances. Only the main 
resonance is characterized by Si = 62, but succes- 
sive application of the generalized Farey operation 
gives distances that are rationally related, that is 
S1/S2 — r/s. Consequently, for all resonances, suc- 
cessive nonlinear mixing through the feedback loop 
can generate only a finite number of new frequen- 
cies, preserving the stability of the corresponding 
response. Also the greater the integers r, s, the 
smaller the stability interval, because more desta- 
bilizing frequencies are added through the feedback 
mechanism. 



4. Discussion 

In Fig. ^ we plot the generalized Farey tree struc- 
ture for med by the predicted thr ee-frequency reso- 
nances [ Cartwright et al., 1999b ]. The frequencies 
of the resonances are obtained by recursive appli- 
cation of the generalized Farey sum starting with 
fi/q and /2 /p. This structure accurately describes all 
the three-frequency resonanc es found in the ex peri- 
ment with phase-locked loops [ Calvo et al., 1999 ] and 
in the numerical simulations with differential equa- 
tions and maps. Successive levels in the tree describe 
the ordering of stability widths in each case. The 
generalized Farey tree structure is thus found to gov- 




Figure 4: (a) Experimental results from an electronic circuit of quasiperiodically forced phase-locked loops 
— a three-frequency devil's staircase. The external frequencies /i and /2 are here fixed at 2100 Hz and 3600 
Hz, equivalent to those in Fig. |l[ We have plotted the third frequency of a three-frequency resonance against 
a control parameter (the DC offset of one of the external forces) for all resonances with plateaux larger than 
a certain size, (b) Shows the hierarchy of three-frequency resonances predicted by the generalized mediant 
starting from the parents /i/4 and ^ijl. At each level in the hierarchy, the daughter resonance formed by the 
mediant between two adjacent parents is seen to be the largest in its interval. 



ern the hierarchy of three-frequency resonances in 
representative dynamical systems with three inter- 
acting frequencies. We conjecture from this that the 
ordering may be universal in three-frequency sys- 
tems. 

From the theorems of Newh ouse, Ruelle, and 
Takens DRuelle & Takens, 197l|; Newhouse et al. 



1978] we expect that this hierarchical structure of 
three-frequency resonances should be relevant to 
the study of torus breakdown and the transition to 
chaos in complex arrays of coupled nonlinear oscilla- 
tors. Such oscillator networks occur in many biolog- 
ical systems, from fireflies and circadian rhythms to 
physiological and neurological systems such as the 
heart and brain. We have investigated the applica- 
tion of these ideas to one such problem in biology: 
that of the mechanism of pitch perception in the au- 
ditory system. We find good agreement between dy- 
namical systems theory and perceptual experiments 
[[Cartwrighte^a/., 1999ap. 
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Figure 5: The power spectrum of the output of the circuit for a DC offset of 0.95 V. Other parameters are as in 
fig. ^. (a) From - 4000 Hz. (b) Detail from 400 - 600 Hz, showing peak at js ~ 518.18 Hz and minor peaks 
separated a distance 6 = efg. 
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